Cosmological Consequences of Exponential Gravity in Palatini Formalism 
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We investigate cosmological consequences of a class of exponential f{R) gravity in the Palatini 
formalism. By using the current largest type la Supernova sample along with determinations of the 
cosmic expansion at intermediary and high-2: we impose tight constraints on the model parameters. 
Differently from other f{R) models, we find solutions of transient acceleration, in which the large- 
scale modification of gravity will drive the Universe to a new decelerated era in the future. We 
also show that a viable cosmological history with the usual matter-dominated era followed by an 
accelerating phase is predicted for some intervals of model parameters. 

PACS numbers: 95.30.Sf, 04.50.Kd, 95.36.-l-x 



I. INTRODUCTION 

Nowadays, there is great interest in modified f{R) 
gravity (see [l[ for recent reviews). Such theories are 
interesting in that they generalize Einstein's general rela- 
tivity and provide insights into the consequences of quan- 
tum corrections to its equations in the high energy regime 
(i? — )■ 00). In cosmology, the interest in these theories 
comes from the fact that they can naturally drive an 
accelerating cosmic expansion without introducing dark 
energy, as happens for instance in the standard ACDM 
cosmology. However, the freedom in the choice of differ- 
ent functional forms of f{R) gives rise to the problem of 
how to constrain the many possible f{R) gravity theo- 
ries. In this regard, much efforts have been developed 
so far, mainly from the theoretical viewpoint [2|. Gen- 
eral principles such as the so-called energy conditions [1| , 
nonlocal causal structure have also been taken into 
account in order to clarify its subtleties. More recently, 
observational constraints from several cosmological data 
sets have been explored for testing the viability of these 
theories 

An important aspect that is worth emphasizing con- 
cerns the two different variational approaches that may 
be followed when one works with /(i?) gravity theories, 
namely, the metric and the Palatini formalisms (see, e.g., 
[l|). In the metric formalism the connections are as- 
sumed to be the Christoffel symbols and the variation 
of the action is taken with respect to the metric, whereas 
in the Palatini variational approach the metric and the 
affine connections are treated as independent fields and 
the variation is taken with respect to both. Because in 
the Palatini approach the connections depend on the par- 
ticular f{R), while in metric formalism the connections 
are defined a priori as the Christoffel symbols, the same 
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f{R) seems to lead to different space-time structures. 

In fact, these approaches are certainly equivalents in 
the context of general relativity (GR), i.e., in the case 
of linear Hilbcrt action; for a general f{R) term in the 
action, they seem to provide completely different theo- 
ries, with very distinct equations of motion. The Pala- 
tini variational approach, for instance, leads to 2nd or- 
der differential field equations, while the resulting field 
equations in the metric approach are 4th order coupled 
differential equations, which presents quite unpleasant 
behavior. These differences also extend to the observa- 
tional aspects. For instance, we note that cosmological 
models based on a power-law functional form in the met- 
ric formulation fail in reproducing the standard matter- 
dominated era followed by an acceleration phase @ (see, 
however, Q), whereas in the Palatini approach, analy- 
sis of a dynamical autonomous systems for the same La- 
grangian density have shown that such theories admit the 
three post inflationary phases of the standard cosmolog 
ical model Q . Although being mathematically more sim- 
ple and successful in passing cosmological tests, we do 
not yet have a clear comprehension of the properties of 
the Palatini formulation of f{R) gravity in other scenar- 
ios, and issues such as solar system experiments Jol , the 
Newtonian limit [lo| and the Cauchy problem [Tlj are 
still contentious (regarding this last issue, see however 
Refs. [13). 

In this paper, we explore cosmological consequences of 
a class of exponential /(i?)-gravity models in the Pala- 
tini formalism. In order to test the observational via- 
bility of these scenarios, we use one of the latest type la 
Supernovae (SNe la) sample, the so-called Union2 compi- 
lation [1^ alon g wi th 11 determinations of the expansion 
rate H{z) [13, H^. In what concerns the past evolu- 
tion of the Universe, we show that for some intervals of 
model parameters a matter-dominated era is followed by 
a late time accelerating phase, differently from some re- 
sults in the metric approach. Another interesting feature 
of this class of models is the possibility of a transient cos- 
mic acceleration, which can lead the Universe to a new 
matter-dominated era in the future. This particular re- 
sult seems to be in agreement with current requirements 
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from String/M theory. 



A. The Gravity Model 



II. PALATINI f{R) COSMOLOGIES 

/(i?)-cosmologies are based on the modified Einstein 
equations of motion derived from the dubbed f{R) grav- 
ity. The action that defines an f{R) gravity is given by 



-gf{R) + S„ 



(1) 



where = SttG, g is the determinant of the metric ten- 
sor and Sm is the standard action for the matter fields. 
Treating the metric and the connection as completely in- 
dependent fields, variation of this action with respect to 
the metric provides the field equations 



I 



while variation with respect to the connection gives 



0. 



(2) 



(3) 



In (m T^i, is the matter energy-momentum tensor which, 
for a perfect-fiuid, is given by T^^ = {pm + Pm.)Ufj.Uu + 
Prndnu, where pm is the energy density, Pm is the fluid 
pressure and is the fluid four-velocity. Here, we adopt 
the notation /' = df /dR and /" = d'^f/dR^. Equation 
give us the connections V^^, which are related with 
the Christoffel symbols {j^j,} of the metric g^jy by 
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Kd^. - 3^..g'"'d,) r . (4) 



Note that in ©, R^j^i, must be calculated in the usual 
way, i.e., in terms of the independent connection F^^, 
given by (|4]) and its derivatives. 

We assume a flat homogeneous and isotropic 
Friedmann-Lemaitre-Robertson- Walker universe whose 
metric is (7^,^ = diag(—l, a^, a^, a^), where a{t) is the cos- 
mological scale factor. The generalized Friedmann equa- 
tion, obtained from ([2|), can be written in terms of the 
redshift parameter z = ag/a — 1 and the density param- 

K^Pmo/(3Fo) as 



H = Hn 



6/' (1 



9 ^ ggn„.o(l + ^) = 

2 /' Rf'-f 



1/2 



(5) 



where p,„o is the matter density today. In terms of these 
quantities, the trace of Eq. ([2]) also provides an impor- 
tant relation: 



Rf - 2/ = -3ij2f^mo(l + zf 



(6) 



There has been an increasing recent interest in models 
of exponential gravity (see, e.g., [l6l - [l9| and references 
therein). Here, we explore a theory of the type: 



f{R) = R- anHl f 1 - g-^/"^") , 



(7) 



in the Palatini formalism. In the above expression, n 
is a free parameter of the theory and a corresponds to 
the strength to which the curvature R scales with the 
Hubble parameter Hq. A functional form of this type 
was originally studied in the metric formalism by Cog- 
nola et al. [ig and Finder [l^l, and some viability con- 
ditions of this model were investigated in Refs. [la, [l9| . 
In the next section, we discuss some cosmological conse- 
quences of the exponential gravity theory given by Eq. 
([7]). To perform our analysis, we impose the positivity 
of the effective gravitational coupling > (to 

avoid anti-gravity and to guarantee that the graviton is 
not a ghost in the sense of a quantum theory), which 
leads to the constraint n < exp (R/aHQ). 



III. OBSERVATIONAL CONSTRAINTS 

Figures 1 and 2 show, respectively, the evolution of 
the Hubble parameter (Eq. [S]) and the predicted distance 
modulus fi{z) = 5 \og[dL (z)] -I- 25 as a function of redshift 
for some best-fit values for n, a and ilm discussed in this 
paper. In the latter expression, = (1 + z) h{z') 
stands for the luminosity distance (in units of mega- 
parsecs and c ~ 1). For the sake of comparison, the 
standard ACDM prediction with f2,„ = 0.27 is also shown 
(thick line). Note that all models seem to be able to re- 
produce fairly well both the H{z) and SNe la measure- 
ments. 



A. Hubble Evolution 

The data points in Fig. 1 arc H{z) determinations 
taken from Ref . [lH • These determinations are based 
on the differential age method that relates the Hub- 
ble parameter directly to the measurable quantity dt/dz 
by H{z) = —z/{l + z) [2^, and can be achieved from 
the recently released sample of old passive galaxies from 
Gemini Deep Deep Survey (GDDS) [2l| and archival 
data 



^ The same data, along with other age estimates of high-2 ob- 
jects, were recently used to reconstruct the shape and redshift 
evolution of the dark energy potential |15| . to place bounds on 
holography-inspired dark energy scenarios [23^ , as well as to im- 
pose constraints on some classes of f{R) models [2^ . 
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FIG. 1: The predicted Hubble evolution H{z) as a function 
of the redshift for the exponential gravity model given by Eq. 
0. The curves correspond to the best-fit values of Qm n 
and a obtained from H{z), Union2 and Union2 + H{z) data 
analyses discussed in the text. For the sake of comparison, 
the standard ACDM model prediction is also shown. The data 
points are the measurements of the H{z) given in Ref. 



FIG. 2: Hubble diagram for 557 SNe la from the Union2 
sample [l^ . The curves correspond to the best-fit values of 
r2,„ n and a arising from statistical analyses involving H{z), 
SNe la and H{z) +SNe la. 



B. SNe la 



In order to impose quantitative constraints from these 
H(z) data on models of exponential gravity, as given by 
Eq. ([7|); we minimize the function 



Nh 

E 

1=1 



(8) 



In the above expresion, Hli^{z\P) is the theoretical Hub- 
ble parameter at redshift Zi, which depends on the com- 
plete set of parameters P = (ri„i,n,a), Hlf^^{z) stands 
for the values of the Hubble parameter given in Ref. [l^l 
and (Ji is the uncertainty for each of the Nh = 12 de- 
terminations oi H{z). In our analysis, we added to this 
H{z) sample a recent estimate of the current value of the 
Hubble parameter, Hq = 72 ± 8 km.s~^.Mpc~^, as given 
by the final results of the HST key project [1^. 

Table I shows the results of our statistical analysis. 
From the above Xh function we construct a likelihood 
function C oc exp(— x|f^/2) and derive the 1,2 and 3a 
intervals for the parameters 57™, n and a. The best-fit 
parameters are the values P that maximize C and the 1, 2 
and 3(T confidence intervals are defined as the sets of cos- 
mological parameters Po- at which the likelihood >C(Pcr) is 
exp(— 1/2), exp(— 4/2) and exp(— 9/2) times smaller than 
the maximum likelihood C{P). From this H{z) analysis, 
the best-fit values found are fi^ = 0.27, n = 1.56 and 
a = 2.56. As expected, due to the current large uncer- 
tainties on the H{z) measurements, we clearly see that 
these data alone do not tightly constrain the values of n 
and a. 



The number and quality of SNe la data available for 
cosmological studies have increased considerably in the 
past few years. One of the most up-to-date SNe la data 
sets has been compiled by AmanuUah et al. [l3|, the 
so-called Union2 sample. This sample is an update of 
the original Union compilation that comprises 557 data 
points including recent large samples from other surveys 
and uses SALT2 for SN la light-curve fitting. 

Similarly to the H{z) test, we estimated the best-fit to 
the set of parameters P by using a statistics, with 



JVsjVe 

XsNe = E 



(9) 



where fj,p{z\P) is the predicted distance modulus given 
above, /^^(z) is the extinction corrected distance modulus 
for a given SNe la at Zi, and <7i is the uncertainty in the 
individual distance moduli. Since we use in our analyses 
the Union2 sample (see [l3[ for details), Nswe = 557. 

1, 2 and Scr intervals for the parameters ^Im, n and 
a arc shown in Table I. For this SNe la analysis the 
best-fit values are fim = 0.33, n = 0.60 and a = 1.48 
with Xmin/'^ — 0-97 stands for the number of degree 
of freedom). Note that the intervals are now consider- 
ably tighter than those obtained from the H{z) analysis 
described above, which reflects the greater constraining 
power of SNe la data when compared with the current 
H{z) sample. 

For completeness, we also performed a joint analysis 
by considering Xt — Xh + x'sNe- The best- fit values for 
this analysis are flm = 0.32, n = 0.56 and a = 1.56 with 
xLn/'^ - l-O- At 2cr level, we found 0.31 < fi,™ < 0.33, 
0.55 <n< 0.92 and 1.53 <a< 1.60 (see also Table I). 
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TABLE I: 1,2 and 3a intervals for f{R) parameters 
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FIG. 3: Deceleration parameter as a function of z for the best- 
fit values of fim n and a presented in Table I. Note that for 
some combinations of parameters the cosmic acceleration is a 
transient phenomenon. The ACDM model, whose predicted 
cosmic acceleration is eternal, is also shown for the sake of 
comparison. 



FIG. 4: Effective equation of state as a function of redshift 
for the exponential gravity theory of Eq. (O in the Palatini 
formalism. As in the previous figures, the curves correspond 
to the best-fit values of ilm n and a obtained from H{z), 
Union2 and Union2 + H{z) data analyses and shown in Table 
I. 



IV. COSMOLOGICAL CONSEQUENCES 

A. Acceleration history 

Refs. [i^ have seriously pointed out a possible con- 
flict between an eternally accelerating universe and our 
best candidate for a consistent quantum theory of grav- 
ity, i.e., String/M theories. The reason is that the only 
known formulation of String theory is in terms of S- 
matriccs, which require infinitely separated, noninteract- 
ing in and out states. As is well known, in the standard 
ACDM scenario the universe will asymptotically become 
a de-Sitter space, which has a cosmological event horizon 
(Ah = / da/a?H{a)) with physics confined to a finite 
region and, therefore, no isolated states. 

In this regard, an interesting feature of the exponential 
f{R) gravity discussed above is the possibility of a tran- 
sient cosmic acceleration with Ah — )- oo. To study this 



phenomenon, let us consider the deceleration parameter 

q{z)=^-^^H\z)^l, (10) 

where a prime denotes differentiation with respect to z 
and H{z) is given by Eq. ([5]). 

Figure 3 shows q{z) as a function of the redshift for 
the three sets of best-fit values obtained in the statisti- 
cal analyses of Sec. III. As can be seen from this figure, 
for some combinations of parameters the Universe was 
decelerated in the past, switched to the current accel- 
erating phase at Zq ~ 1 and will eventually decelerate 
again at some Zd < —1. For these sets of parameters, 
it is possible to show that Ah — > oo, thereby alleviating 
the potential theoretical and observational confiict dis- 
cussed above. It is worth emphasizing that this kind of 
dynamic behavior is not found in most of the /(i?) cos- 
mologies discussed in the literature [l^] , being essentially 
a feature of the so-called thawing and hybrid quintessence 
potentials [l^ , some classes of coupled quintessence mod- 
els [2^ and brane- world scenarios . 
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B. Effective equation of state 

In Rcf. 6], it was shown that f{R) derived cosmolo- 
gies in the metric formahsm cannot produce a standard 
matter-dominated era foUowed by an accelerating expan- 
sion (we refer the reader to for a different conclusion). 
To verify if this undesirable behavior happens in the Pala- 
tini f{R) gravity discussed in this paper, we derive the 
effective equation of state (EoS) 

Weff = -l + ^^^H'iz) (11) 

as a function of the redshift. 

In Figure 4, we show the effective EoS as a function of 
z for the best-fit values discussed in the previous section. 
Note that, for some of these combinations of parameters 
(basically those derived from SNe la data), the universe 
went through a past matter-dominated phase {w = 0) be- 
fore switching to a late time accelerating phase {w < 0). 
In particular, we note that for the best-fit values derived 
from the SNe la plus H{z) joint analysis, there seems 
to be evidence for a slowing down of the cosmic accel- 
eration today, which is somewhat in agreement with the 
results of Ref. [sif for some dark energy parameteriza- 
tions. From the results shown in Fig. 5, we clearly see 
that the arguments of Ref. Q about the behavior of w^j / 
in the metric approach seems not to apply to the Pala- 
tini formalism, at least for the exponential f{R) gravity 
theory studied here and the interval of parameters 
n and a given by our statistical analyses. 



V. CONCLUDING REMARKS 

Cosmological models based on /(i?)-gravity may ex- 
hibit a natural acceleration mechanism without intro- 
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ducing a dark energy component. In this paper, we have 
investigated cosmological consequences of a class of expo- 
nential /(i?)-gravity in the Palatini formalism, as given 
by Eq. ([7]). We have performed consistency checks and 
tested the observational viability of these scenarios by 
using one of the latest SNe la data, the so-called Union2 
sample with 557 data points and 11 measurements of 
the expansion rate H{z) at intermediary and high-z. We 
have found a good agreement between these observations 
and the theoretical predictions of the model, with the re- 
duced Xmin/'^ — 1 the three tests performed. 

Differently from the dynamical behavior of other f{R) 
scenarios discussed in the literature (either in metric or 
Palatini formalisms), we have found solutions of tran- 
sient cosmic acceleration in which the large-scale modifi- 
cation of gravity will drive the Universe to a new matter- 
dominated era in the future. As mentioned earlier, this 
kind of solution is in full agreement with theoretical re- 
queriments from String/M theories, as first pointed out 
in Ref. 

Finally, we have also shown that, differently from the 
results of Ref. for power-law f{R) gravity in the met- 
ric formalism, exponential f{R) models corresponding to 
the best-fit solutions from SNe la and SNe la -I- H{z) 
minimization have the usual matter-dominated phase 
followed by a late time cosmic acceleration (see also Q 
for a discussion). 
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